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Abstract 

Fluctuations of the comoving curvature perturbation with wavelengths larger than the horizon length 
are governed by a Langevin equation whose stochastic noise arise from the quantum fluctuations that 
are assumed to become classical at horizon crossing. The infrared part of the curvature perturbation 
performs a random walk under the action of the stochastic noise and, at the same time, it suffers a 
classical force caused by its self-interaction. By a path-interal approach and, alternatively, by the 
standard procedure in random walk analysis of adiabatic elimination of fast variables, we derive 
the corresponding Kramers-Moyal equation which describes how the probability distribution of the 
comoving curvature perturbation at a given spatial point evolves in time and is a generalization of the 
Fokker-Planck equation. This approach offers an alternative way to study the late time behaviour 
of the correlators of the curvature perturbation from infrared effects. 
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1 Introduction 



Cosmological inflation [T] has become the dominant paradigm within which one can attempt to un- 
derstand the initial conditions for Cosmic Microwave Background (CMB) anisotropics and structure 
formation. In the inflationary picture, the primordial cosmological perturbations are created from 
quantum fluctuations which are "redshifted" out of the horizon during an early period of acceler- 
ated expansion. Once outside the horizon, they remain "frozen" until the horizon grows during a 
later matter- or radiation-dominated era. After falling back inside the horizon they are communi- 
cated to the primordial plasma and hence are directly observable as temperature anisotropics in the 
CMB. These anisotropics have been recently mapped with spectacular accuracy by the Wilkinson 
Microwave Anisotropy Probe (WMAP) [2]. 

These CMB observations show that the cosmological perturbations are very small, of order 10^^ 
compared to the homogeneous background. Therefore, one might think that first-order perturbation 
theory will be adequate for all comparison with observations. However, that may not be the case; 
the Planck satellite [3] and its successors may be sensitive to non-Gaussianity in the cosmological 
perturbations. Such non-Gaussianities are sourced by self-interactions in the early universe, and 
become visible at the level of second- or higher-order perturbation theory [4J. This possibility is 
of considerable interest in its own right and is presently being vigorously explored. However, there 
are other equally compelling reasons to go beyond linear theory. Self-interactions of any scalar 
field during the inflationary stage, and more interestingly of the comoving curvature perturbation 
C, give rise also to corrections in all correlators, including the observationally interesting cases 
of the two- and the three-point correlation functions. Such corrections are associated with so- 
called "loops," in which virtual particles with arbitrary momentum are emitted and re-absorbed 
by the fields which participate in the correlation function. Loop corrections to the correlators 
of cosmological perturbations during inflation have been extensively studied in the last few years 
O El m El m [ini E] ■ in particular, it has been known for a long time that there exists a breakdown in 
the perturbative expansion due to infra-red (IR) divergences for the case of an interacting test scalar 
field in a pure de Sitter stage [121 US] • Recently the IR issue was anlyzed in the more interesting (but 
also more challenging) case of real metric perturbations, taking also into account the fluctuations of 
the background, and it was shown that in this case there appears to be a perturbative breakdown 
globally in the "large box" on a time scale, t ~ RS, given by the de Sitter radius and entropy 

, which indicates some interesting connections to black hole physics and the information paradox 

[MKllKilKlT]. 

The semi-classical one-loop IR corrections depend on the (logarthmic value of) infra-red comoving 
momentum cut-off L~^. and there has been quite a debate in the literature on the interpretation of 
these IR corrections [IlIIlEOlEIlEaiMllMlESlESlEHlEHlI^ 

One option will be to restrict ourselves to the computation of cosmological observables within 
our local patch. This means that all computations should be done within a comoving box of size 
L whose present size is not much larger than the present horizon Hq^, in order that the answer 
depend only on local properties of the universe at our position. In such a case, IR corrections are 
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of the order of ln(/cL) = 0(1) if we observe the properties of the cosmological correlators on a given 
observationally interesting comoving wavenumber k. 

The other option is to set L to be some comoving length scale which left the horizon many 
e-folds before the observable universe. In particular, the smallest possible value of is aiH, 
where Oi is the value of the scale factor at the beginning of inflation. Since the wavelength k^^ 
goes out of the horizon when the scale factor equals a^, = k/H, in such a case ln(A;L) = ln(afc/a,) 
is proportional to the total number of e-folds from the beginning of inflation to the time when the 
mode k exits the horizon. Nevertheless, these large IR corrections appear only in quantities which 
are not directly observable [191 EI] today. Indeed, one should be suspicious that large effects from 
IR modes are somehow unphysical for a present observer, because on general grounds we expect 
the structure within our observable patch of the universe to depend only on its local properties. 
Although their appearance cannot be avoided, they do not themselves have any particular interest 
unless one wishes to ask more general questions, for instance, what is the probability to find an 
inflaton field homogeneous enough to lead to the correct temperature anisotropy in the CMB within 
patches like ours? Or, what is the probability that the power spectrum of C has the value we 
measure? If one insists in adopting a superlarge box, it should be kept in mind therefore that the 
location of our box may be untypical and that the large IR effects appear as a measure of how 
likely it is that the correlators averaged over the superlarge box coincide with the correlators in our 
observable universe. One may also be interested in understanding the perturbative breakdown in 
the large box, because it teaches us something fundamental about de Sitter space-times, which can 
be valuable for understanding also other problems. In any case, there appears to be a growth of 
perturbations in de Sitter, and we are interested in properly understanding its physics. 

To deal with this problem, but also for other less ambitious reasons, one can try to use the 
approach of stochatic inflation to study the behaviour of the IR part of the comoving curvature 
perturbation ^. The stochastic approach embodies the idea that the IR part of the field may be 
considered as a classical space-dependent stochastic field satisfying a local Langevin-like equation 
|33j . The stochastic noise term arises from the quantum fluctuations which become classical at 
horizon crossing and then contribute to the background. The IR part of the field C performs a 
random walk under the action of this stochastic noise and, at the same time, it suffers a classical 
force caused by its self-interaction. Physically, the stochastic force is caused by the ultraviolet 
modes that are instantaneously redshifting across the horizon at time t. It is Markovian because 
each increment is suitably chosen to be uncorrelated with the previous one. As the stochastic 
approach seems to reproduce the terms with the largest power of the IR logarithms at each order 
in the coupling constant (leading log approximation) [HU |35] , it is therefore useful to study the late 
time contributions of the correlation functions of C and to assess whether these can be computed by 
a classical approximation. 

The interaction terms of the comoving curvature perturbation ( are purely derivative type inter- 
actions, which has for some time hindered a progress in solving the stochastic equations for In the 
present paper, we will attempt to make some modest progress on that front. We believe that this can 
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shed some additional light on the IR issues, although a very simple different semi-classical approach 
to study IR loop corrections to the comoving curvature perturbation and the tensor perturbations 
has been developed in Ref. [H] (see also [271 ES] for further discussions of this approach) . 

For small internal momenta, one can derive the interaction terms for the various degrees of 
freedom in a specific model of inflation, and then use the conditions of the theorem of Ref. [5] 
to decide whether these interactions can lead to anomalous late time contributions. For single 
field inflation (possibly together with free massless scalar fields cj's coupled to it) it was shown in 
Ref. [5] that the interactions do obey the conditions, and therefore do not give anomalous late 
time contributions to all orders. This can be compared with [36], where it was shown that C is 
indeed conserved after horizon exit as a consequence of the fact that the system possesses a time 
translational invariance in e-folding time when time derivatives dominate over spatial ones. The 
stochastic approach (which in any case has nothing to say about the ultraviolet loops) may offer 
an alternative way to study the late time behaviour of the correlators of ^ and to show that the 
resummed IR loop corrections from these modes that crossed the horizon arbitrarily in the past do 
not give rise to any additional time dependence of the two-point function of the comoving curvature 
perturbation, which therefore stays constant once the mode k is outside of the horizon. 

The goal of this paper is to obtain the equation governing the probability P{C, t) that, at a given 
point, the comoving curvature perturbation acquires a given value (" at a given time t, in the presence 
of self-interactions. Some aspects of this subject were investigated in a seminal paper by Bond and 
Salopek [37] and more recently in Refs. [381 EH HOI HIl HH SS]- As we shall see, the equation for 
P(C, t) is not a diffusive-like equation, that is a Fokker-Planck equation, but a generalization of it, 
the so-called Kramers-Moyal equation. It contains an infinite number of derivative with respect to 
C (the Fokker-Planck equation contains up to two derivatives). This is an unavoidable consequence 
of the fact that self-interactions of the comoving curvature perturbation are present. A key point 
we would like to anticipate about the Kramers-Moyal equation is that it will prove to be convenient 
to derive it starting from the probability P{C,C,t), i.e. in terms of two variables, C and C. The 
stochastic motion will then be characterized by the fact that C decays rapidly, within a few Hubble 
times, to an equilibrium value due to the presence of a nonvanishing potential V{C). As we will 
be only interested in times larger than a few Hubble times, the fast variable C will be integrated 
out, leaving us with an effective equation for the probability of the slow variable C. This adiabatic 
procedure is similar to what is done in field theory where heavy degrees of freedom are integrated 
out resulting in an effective Lagrangian describing the dynamics of the light fields. The fact that C 
rapidly settles to its equilibrium value allows the effective diffusion constant of the C field to reach 
a constant (in time) value, thus avoiding an anomalous growth of the cumulants of C at large times. 

The paper is organized as follows. In section 2 we derive the generic, but exact Kramers-Moyal 
equation from the a path-integral formulation. In section 3 we will discuss an alternative derivation 
of the Kramers-Moyal equation through the adiabatic elimination of fast variables in the case in 
which the self-interaction of the curvature perturbation is dominated by of a cubic time derivative 
nature. Section 4 deals with the case in which the curvature perturbation is coupled to an extra 
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scalar degree of freedom. All these examples are studied under the assumption of neglecting the 
time evolution of the slow-roll parameters. In section section 5 we relax this assumption and study 
a particular case of slow-roll single field inflation. Finally, section 6 offers our conclusions. 



2 The Kramers- Moyal equation for the curvature 
perturbation 

The scope of this section is to describe a generic method to obtain the exact equation governing 
the probability of finding a given value of the comoving curvature perturbation ^ at a given time t. 
The method is based on the path-integral formulation whose technical details can be found in Refs. 



Before analyzing in detail the case of the comoving curvature perturbation, we will first discuss 
the more standard example of a self-interacting scalar field $ with potential V{^) in de Sitter space, 
i.e. we assume that the Hubble rate H is constant in time. Its field equation is given by 



^ + 3H^- irV^^ + V'(<^) = 0. (1) 



The stochastic approach is based on the fact that the IR part of the field may be considered as a 
classical field satisfying a Langevin-like equation [33]. To find the equation for the infrared part of 
the field 

^ e {Ha - k) [e^-^MtW + e-^'^-"./'j;(t)aj,J , (2) 
one can follow the rules rigously obtained Ref. [35] applied to Eq. ([1]): 

• At each order in the field, retain only the term with the smallest number of derivatives; 

• for the linear terms in the field, each time derivative has a stochastic source substracted. 
This set of rules gives 

4) + - ^V^cf) ^ 3H<j) ^3h(^^- r?^) , V'{(j)) ^ ^'(0) , (3) 

where 

^ kSDik- aH) Y''-^Mt)a^ + e'^''-^.^^ (t)aj,J , (4) 

is a Gaussian stochastic noise. Its two-point correlator computed over many ensemble realizations 
is 

(7?^(x,t)?7<^(x,t')) = -^^D{t - t') . (5) 
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The behaviour of the IR part of the field therefore is described by a Langevin-hke equation 

<^(x,t)=r?^(x,t)-^. (6) 

The IR part of the field (p performs a random walk under the action of the stochastic noise and, 
at the same time, it suffers a classical force caused by its self-interaction. Physically, the stochastic 
force is caused by the ultraviolet modes that are instantaneously redshifting across the horizon at 
time t. It is Markovian because each increment is uncorrelated with the previous one. This is due to 
the choice of the step function in momentum space to select the IR modes. Any other choice of the 
smoothing function would lead to a non-Markovian noise and to memory effects in the stochastic 
motion which will require a proper treatment [44]. 

In order to obtain the Fokker-Planck equation describing the probability t) that the IR 
part of the field acquires a given value (/) at a given time t, we exploit the path-integral formulation. 
We consider an ensemble of trajectories all starting at t = from the initial position <?!<(x, 0) = </>o 
and we follow them for a time t. We discretize the interval [0, t] in steps At = e, so t^ = ke with 
k = 1, . . . n, and tn = t. A trajectory is then defined by the collection of values {(/>!,...,(/)„}, such 
that (p{tk) = 0fc (we do not indicate anymore the spatial variable from now on). The probability 
density in the space of trajectories is 

W{cj)o; (/>!,..., <Pn; tn) = {SDiHh) - <Pl) ■ ■ ■ <5D(</'(in) " </'n)> • (7) 

Our basic object will be 

/oo fOO 
dcl)i... d(i)n~lW{(l)o;(pi,...,(pn--l,<Pn]tn)- (8) 

-oo J —oo 

Using the integral representation of the Dirac delta 

Sd{x) = I ^e-^^^ (9) 

we write Eq. ([7]) as 

J ZTT zvr 

We must therefore compute 

= (e-^^-i^'<^(*^)). (11) 

This is a well-known object both in quantum field theory and in statistical mechanics, since it is 
the generating functional of the connected Green's functions, see e.g. [46j. To a field theorist this is 
even more clear if we define the "current" J from —iX = eJ, and we use a continuous notation, so 
that 

^ (gi/dtj(t)0(t)^ ^ (12) 
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Therefore 



n n 



1 

-iXi{(j)i)c - {(pi4>j)c H ^ AjAjAfc {(j)i(pj(l)k)c 

H — ^ XiXjXkXi {4'i4>j(t)k4>i)c + • • • , (13) 



where the sum over i, j, ... is understood. This gives 



poo 




poo 




/ ' 


lct>^.. 






/ — oo 




J — oo 





(14) 

To get the equation that governs the evolution of the probability P{(pn, tn), we take the time deriva- 
tive of Eq. ([IID 

J-oo J 27r 27r 



dP_ ^ _.d(Mc 

dtn dtn 

- ^-^-^^^rdct>,...d<l.n-J^-^...'^-^Xkk^^^^^^^ 

2 dtn J-oo J 27r 27r 

+ •••• (15) 

Using the replacement Aj exp{iAi(/)i} = — iS^^ exp{iAi(/)i} (the index z is not summed over), inside 
the integrals we can replace Aj — —id^^. Since we integrate over dcpi, . . . d(j)n-i, but not over dcpn, if 
the generic index n the term 9^-, when integrated over dcpi, is a total derivative and gives zero, 
because at the boundaries (pi = ±00 the integrand vanishes exponentially, and the only contribution 
comes from i = n. The equation (fT5]) reduces therefore to (setting tn = t and 0„ = (/>) 



ap(0, t) _ (-ir g/ip(t) o^p(0,t) 



dt »! 5t 



where 



/ip(t) = (0?'(t)), (17) 

is the p-th connected cumulant. Equation (|16p is called the Kramers-Moyal (KM) equation or "the 
stochastic equation" , and is well known in the theory of stochastic processes [Ml [47] . Its solution 
requires the exact knowledge of the cumulants and in this sense the KM equation does not contain 
more information than the generating functional. Furthermore, Pawula's theorem [50] states that if 
we assume fip = for some even p, then we are actually assuming fip to be zero for all p > 3. Thus, 
it is logically inconsistent to retain more than two terms in the KM expansion unless all of the terms 
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are retained. Nevertheless inconsistencies may lead to tiny numerical errors and one usually reduces 
the KM equation (|16p to a Fokker-Planck equation 



dt 



J_d_ 
m'd4> 



P{<l>,t) 



(18) 



where one has exploited the fact that at leading order in the coupling constant (s) the connected 
two-point correlator is dominated by the Gaussian and Markovian stochastic noise. At small times 
the diffusive term dominates 



and the scalar field performs a diffusive random walk such that {(f)'^) ~ H^t. At larger times, the 
friction term intervenes and the equation might have a non-perturbative stationary late time solution 
induced by the fact that the scalar force eventually balances the tendency of inflationary particle 
production to force the scalar up its potential [H] 



To estimate the time scale at which the stationary solution is reached one has to specify the potential. 



the field $ is the inflaton, the eternal inflation regime is achieved when the potential is flat enough 
that quantum corrections always dominate over the classical motion. In such a case, the diffusive 
term dominates and the variance of the field grows linearly in time. 

The derivation of the equation governing the probability P{C, t) through the path-integral formu- 
lation follows basically the same steps described for the self-interacting scalar field in de Sitter space. 
In order to find a simple solvable toy model, we find it convenient to adopt the effective field theory 
of inflation [19] which makes it possible to explore in full generality all the possible self-interactions 
of the inflaton, and in particular to clearly see that the inflaton can have large self-interactions 
without spoiling the background quasi de Sitter solution. It was shown that the self-interactions 
can be parametrically larger than the one mediated by gravity and therefore we can in some cases 
concentrate on those without worrying about metric fluctuations. Note, however, that in doing this, 
we are exactly neglecting the interactions that leads to large IR effects. The interactions that leads 
to large IR effects are subdominant in the slow-roll expansion, and vanishes in the decoupling limit 
of the effective approach. So, while we don't expect to see any large IR effects in this limit, it is still 
useful for providing a simple solvable toy model for testing the stochastic approach. 

The effective Lagrangian is a function of the fleld vr which represents the Goldstone boson of 
time translations which are spontaneously broken during inflation. It is related to the (" metric 
fluctuation at linear level by the simple relation = —Htt. The general derivation of the Lagrangian 
for vr can be found in Ref. 09] and we do not present it here. It suffice to say that, for instance, 
at the cubic level there are in general two interactions of comparable strength, tt'^ and 7r(V7r)^, each 




(19) 




(20) 



For instance, for a quartic potential V{(l)) = Xc/)'^, the variance of the field is given by {(/P) ~ H'^/y/X. 
This means that it takes a time t ~ /\/\ for the stationary regime to be reached. Notice that if 
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one accompanied by quartic interactions due to symmetry reasons. Upon tuning, one can make the 
term containing vr^ parametrically larger (smaller) than the one containing 7r(V7r)^. Following the 
notation of [49], the action for reduces to the form (see also |51j ) 



.g|,M?,(c=-|(vc)^ 



PI 



4 M| -3 
_l lA'^ _|_ 



(21) 

where Cg is the sound speed, H = —eH^ defines the slow-roll parameter e, M3 is a mass term and the 
dots stand for the terms which are proportional to the slow-roll parameters and to the higher-order 
terms in perturbation theory. The IR part of the field can be defined as usual 



C(x,t) 



(2vr) 



3 . (Fa - k) e*'^-"Ck(t)ak + e-'''-"Ck(*)4 



It is subject to a stochastic noise 



d^k 



kSoik- aH) (^^■''C\.{t)a\, + e-^'^'^CklOa; 



k-x A* / 



(22) 



(23) 



(27r)3 SKVV-K , - skv^;-kj 

whose two-point correlator computed over many ensemble realizations is computed using linear 
perturbation theory and is defined as 



D. 



where 



SD{t-t') 



(24) 



(25) 



The corresponding motion is stochastic and the probability that the curvature perturbation C ac- 
quires the value after n-steps at time tn is given by 



00 fOO 

dCi.. 



/ 27r 2tt 



(26) 

Taking the derivative with respect to t„ and setting at the end tn = t and Cn = C we obtain the 
exact Kramers-Moyal equation for the comoving curvature perturbation 



dp{c,t) ^j^{-iy dfip{t) dppict) 



dt 



p=2 



pi dt dQP 



(27) 



where 



^^p{t) = {em 



(28) 



9 



is the p-th connected cumulant of the perturbation ^. The sum over the index p is restricted to 
p > 2, the term proportional to one derivative with respect to C, is not present. This is because 
we require that (, being a perturbation, has an expectation value which vanishes at any order in 
perturbation theory. If we start with the initial condition (C)c = at t = 0, then equation ([27]) 
insures that it remains zero at all times 

|(C)c = I dCC§-^PiCt) = o, (29) 

as one can see by plugging in expression (j27|) and integrating by parts. The KM equation (p7|) for the 
curvature perturbation does not contain any more information than the full generating functional of 
In this sense, there is complete analogy with what happens in field theory where the knowledge of 
the full dynamics of a system requires the knowledge of all the n-point connected correlators. The 
KM equation describes the time-flow of the functional in the same way the renormalization group 
equations describe the flow of the generating functional and the n-point correlators with energy. As 
in field theory where one proceeds performing a perturbation expansion in some coupling constants, 
one can try to perturbatively simplify the KM equation to solve for the probability. Even if Pawula's 
theorem [50] implies that throwing away terms in Eq. (|27p with a number of derivatives with respect 
to C higher than two to obtain a diffusive Fokker-Planck equation is logically inconsistent, it might 
lead to tiny errors [52] , especially if one considers the case of the evolution during the last 60 e- folds 
of inflation. For a cubic interaction one could make an expansion in the coupling constant(s) keeping 
only the cubic term p = 3 and approximating the two-point cumulant to its tree-level value deduced 
from (f24]) . 

(C\ = /'d*'^. (30) 

Alternatively, one can stop at the second derivative with respect to C and obtain a Fokker-Planck 
equation. Again, as in field theory the perturbative renormalization group evolution of the couplings 
constants is able to catch and resum the leading logs, the KM equation may be able to provide non- 
perturbative solutions, like the stationary one in Eq. ()20p . As far as the initial condition is concerned, 
Eq. (p7|) should be solved with the assumption that at t = 

P{C,0) = 5d{C)- (31) 

Let us pause for a moment and discuss the interpretation of such initial condition \37\ [21] . If we 
consider only the physical observables within a comoving region coinciding with our present horizon 
volume, where several observations indicate that this patch is highly homogeneous and isotropic, 
we have to select those initial conditions which respect these observations, and therefore select the 
initial time t = t^Q to be the time when one is left with the last 60 e-folds or so of inflation. If, 
on the other side, we insist in adopting a superlarge box containing many smaller patches of size 
similar to our observable universe, as commented in the introduction, the initial time t = should 
be taken to be the beginning of inflation. If so, it should be kept in mind that the location of our 
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box may be untypical and one should quantify how likely it is that the correlators averaged over 
the superlarge box coincide with the correlators in our observable universe. One should split the 
computation of the probability in two steps. First, the distribution of the values of the curvature 
perturbation inside the superlarge box should be estimated. This will allow one to compute the 
probability P^Ceo I ^6o) that the starting initial condition for the background field is given by Ceo in 
our local observable patch at time igo when there remain 60 or so e-folds till the end of inflation 
in that patch. The distribution of the comoving curvature perturbation in the superlage box, due 
to the IR divergences, will be far from being Gaussian. However, we do expect that in our local 
patch initial conditions are relatively uniform. From the point of view of our local universe the IR 
properties of the superlarge box only enter to provide the probability -P(C6o|^6o) that the value Ceo 
is achieved in our patch. At later times the field value at a given point has a distribution given by 



Furthermore, if the comoving curvature is homogeneously distributed in our local patch, then 
P{C I i6o;i60,C6o) = '^d(C-C6o) with Cgo - 0. The probability P(C,i) is generically highly non- 
Gaussian, because the initial conditions are not Gaussian distributed as large IR corrections are 
present in the superlarge box. 



3 An alternative derivation of the Kramers- Moyal equa- 
tion: the adiabatic elimination of fast variables 



The scope of this section is to offer an alternative and perturbative derivation of the KM equation for 
the curvature perturbation. We will restrict ourselves to the case in which large non-Gaussianities 
are present, i.e. we will disregard the self-interaction terms in the Lagrangian proportional to the 
slow-roll parameters. The same assumption allows us to utilize the linear relation between C and vr. 
Again, we note that we don't expect any large IR effects in this limit, since the interactions leading 
to large IR effects are slow-roll suppressed, but it may nevertheless be a useful limit for testing the 
methods. 

We are interested in the behaviour of the long wavelength part of the comoving curvature per- 
turbation. Following the procedure described in the previous section, we consider the IR part of C 
as a classical space-dependent stochastic field which satisfies a local Langevin-like equation. Fur- 
thermore, we suppose that at the cubic level the interaction term is dominated by the term, i.e. 
we suppose that M| » eff^Mpj. Following the rules dictated by [35], the IR part of the comoving 
curvature perturbation satisfies the equation 




(32) 



where r/^ is the stochastic Markovian noise whose two-point correlator i 




(33) 
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(7?c(x,t)r/c(x,t')> = ^^J^]^^^D{t-t') . (34) 

To determine the equation governing the probabihty, it is useful to rewrite Eq. ([33]) doubhng the 
degrees of freedom [U] 



C = v^, (35) 
+ 3Hv^ + 6cl-^t_vl = 3i/77c ■ (36) 



Fohowing the steps described in the previous section, we may write the KM equation for the prob- 
ability P{(,Vi^,t) as 



where 



Viv,) = '-Hvl + 2cl-S7T-l (38) 



PI 

Notice that in the first and second terms of the right-hand side of Eq. (j37p we might replace 
by (u^ — {v(^)) and V'{v(^) by {V'{C) — (l^'(w^))) to account for the fact that the comoving curvature 
perturbation ( and must have vanishing expectation value at all times at all orders. We will come 
back to this point later. The initial condition of Eq. (I37p is given by 



P{C,v^,0) = 5d{0Sd{vc) • (39) 

In order to solve Eq. ()37p we may proceed using the method of the adiabatic elimination of the 
fast variables [54] • The idea is that the variable will decay rapidly to an equilibrium value due 
to the presence of a nonvanishing potential l^(f^). Inspecting the energy scales in the potential, it 
is easy to get convinced that the equilibrium situation should be reached after one Hubble time, 
icq ~ H^^. As we are only interested in times larger than the decay time teq of the fast variable 
V(^, the process described by Eqs. (j35|) and ([36]) is then mainly described by the motion of the slow 
variable Hence we may say that the slow variable slaves the fast variable [55]. The slow variable 
is then the important or relevant variable, since the fast variable becomes irrelevant for the above 
time scale because it can be expressed by the slow variable. Integrating out the fast variable Vt^ 
will allow us to get an effective equation for the probability of the slow variable ^ with a procedure 
similar to what is done in field theory where heavy degrees of freedom are integrated out resulting 
in an effective Lagrangian describing the dynamics of the light fields. 
We first look at the operator 
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We introduce the following operators 




2 dv^ y 2D^ ' 



which, for natural boundary conditions are the adjoints of each other, i.e., a = . The operator 
L^^ can be recast in the form 

= -exp {-V/D(;) a^a exp {V/D^) . (42) 
If we denote by i^niv^) the eigenfunctions of the Hermitian operator a 

a^aipn = Kipn, (43) 

then all A„ must be real and non-negative. Furthermore = exp {—V/D(^) ipn are the eigenfunctions 
of the operator L^^ with eigenvalues — A„ < 0. For n = we have the stationary solution induced 
by the zero mode 

*o(^^c) = i-^V/D^) , Ao = , (44) 

where A/q is the normalization factor. The eigenfunctions = V^oi^o and the eigenfunctions 
= ipn/V^ = ^n/^o of the adjoint operator lI^ form an orthonormal basis satisfying 

'^i{V(;)'$miVi;)dV(; = 6^n, (45) 

J — oo 

and 

oo 

Y^^I{v^)^^{v'^) = 6d{vc-v[) . (46) 

n=0 

Notice in particular that ^o('^c) ~ '^^^ expand the distribution function P((,Vi^,t) into the 

complete set of the operator Ly^ 

oo 

P{C,VC,t) = PmiCmmivc). (47) 

m=0 

We insert this expansion into the KM equation (|37p multiplying the resulting equation by and 
then integrating over the variable f^, thus obtaining 

g-^+Xn) Pn=Y. Ln,mPm (48) 
^ m=0 
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with 



Ln,m= '^i{V(;)L(;'^m{v<;)dVt;, (49) 

and 

^C = -^^C- (50) 

The infinite set of equations ()48p is exact and, as expected, it contains an infinite number of deriva- 
tives with respect to C once an iterative procedure is adopted to solve for the Pm^s. At times 
t S> teq ~ H~^, we may use the following approximation procedure. Knowing that the eigenvalues 
Xn are positive for all n > 1, we may neglect in (jl8|) the time derivative in the equations for n larger 
than zero, that is 

g CXI OO 

^-Pq = Lo^m P-m = + -^0,m Pm , 

m=0 m=l 



OO 



dt 

Pn — ^ ^ I-'n,m Pm 

m=0 

OO 

= A„ ^ LnQ Pq + Ln^m Pm 

OO 

— Ki^ ^n,0 Pq + ^n^^^ ^ Ln,m LmflPo + • • • j (51) 

m=l 

where the last passage is possibile because the Pm for m > 1 decays in time faster than Pq for which 
Ao = 0. We thus obtain that the equation governing the probability 

/OO °° fOO 

-OO n=0-^-°° 

^ POO 



n=0 ■ 

PoiCt), (52) 



is the KM equation 
d_ 
di 



i-P{C,t) = LKuPiCt) 



E '^"^ ^^'"^ \ "^"'^ ~'~ '^"^ E ^n,mLmfi H 1 PiC^t) 

n=l \ m=l J 
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which contains an infinite number of derivatives with respect to C. The first two coefficients of the 
KM expansion are 



2 9iJ2 
and 



oo 



I ^eff oo 

I ,T, (y^) dvA{ I Mv^) dv^ ] , (54) 



n=l ^-^-"^ 



_ 5f = ^ U v^^n{vc) dvA ^ n ^iiv^) ^m{v^) dvA U ^lv^^o{v() dvA 

(55) 

Notice that the term proportional to the first derivative with respect to ^ coming from Lo,o vanishes 
which is justified by requiring that at ah orders in perturbation theory the expectation value of C 
must vanish. Indeed, we should replace the operator = v^d/dC with (f^ — {v())d/dC- This implies 
that Lo,o should be replaced by 

(POO POO POO oo \ f) 

/ dv^^lv^^o- dC dv^^PniC'X'fn] ^ 

J —OO J —CO J —oo n=0 / 

/ f'oo f'oo poo f'oo roQ 

/ dv^^lv^^o- dC'PoiC') dv^^lv^^o- dC'PiiO dv^^lv^^i + -- 

\J —oo J —OD J —oo J —OD J —OO 



= 0, 

where we have used the fact that ^'J = 1 and PQ{(,t) = P{C,t) is normalized to unity, and that 
the P„'s with n > 1, according to (jST]) . are expressed in terms of derivatives of Pq with respect to C 
and therefore their total integrals vanish. Due to the Pawula's theorem [50], we are not allowed to 
disregard the terms with a number of derivatives with respect to ^ higher or equal to three in the 
KM equation. Nevertheless, if we are only interested in the time behaviour of the variance (C^) 

|(a=/<^CC^|p(C,.). (57) 

we can restrict ourselves to the term with two derivatives with respect to C- This gives 

(0 = ^^- (58) 

The linear dependence originates by our initial assumption that the slow-roll parameters are tiny so 
that their time evolution is negligible; otherwise we should use as time variable not t but = In a 
[33] and get 

pN tt2 

iC) = / dN' ^ . (59) 



15 



The coefficient D^^ in (I54p gets contribution only at second order in the perturbative parameter 
A = (2c^M|/e-ff^Mpj) in front of the cubic term in the potential (I38p . Indeed, the 0{X) correction 
to the eigenvalues A„ vanishes, 

5A„ = A(^'Wt|^;3|xI/(0)) = 0, (60) 

where "^n^ are the eigenfunction of the linear problem (see below); the correction to the eigenfunc- 
tions are given by 



^^>W = ^rW + 0(A) \^/n{n-l){n-2)^ 



+ 



9^3/2^(0) ^ - 9(n + l)3/2<°ji - V(r^ + l)(n + 2)(n + 3)<°j3l . (61) 



Inspecting the expression (|54p , we can first insert the expression ()6ip in the first integral in the right 
hand side while anywhere else we pick up the eigenfunctions at the zero-th order. This selects an 
odd n from the second integral since is odd and ^'q'^^ is even. However, from (j6ip we read that the 
corrections to the the eigenfunctions with n odd receive corrections only from even functions "^^^ 
and therefore this correction vanishes. We now perturb the second integral in the expression (j54p . 
Being proportional to ^^^^ (see below) this selects n = 1 in the first integral; again the corrections 
in "^^^ are all even and this makes the second integral vanishing. The same happens if we perturb ^'o 
in the second integral as it gets corrections only from odd functions. We conclude that the effective 
diffusion coefficient receives corrections only at second order in the coupling constants. On the other 
hand it is easy to see that S^^ receives contributions which are 0{X). This is what we expect from 
computing diagramatically the one-loop correction to the connected two-point cumulants (C^) and 
{(^^) . The corrections to the diffusion coefficient are time-independent. This is intimately linked 
to the fact that the reduction to the KM equation (j53p has been made possible because the fast 
variable reaches rapidly a stationary solution and can be integrated out. This conclusion is in 
fact valid at any order in perturbation theory when the dominant contribution to the non-linearities 
at any given order p is dominated by the time-derivative interaction = CJ'. In such a case, the 
potential for the fast variable is of the form 

y(-c) = lHvl + ±'-^v^, (62) 

m=l 

and the elimination of the fast variable applies. This allows us to conclude that comoving curvature 
perturbation performs a diffusive random motion where (C^) ~ t (neglecting the time variation of 
the slow-roll parameters) even in the presence of time-derivative interactions: after a time scale of 
the order of the inverse Hubble time the diffusive coefficient saturates to a constant. In the 
opposite case in which the ^ (VC)^ cubic term dominates over the term, one expects an equation 
for the probability which is very similar to the Kardar-Parisi-Zhang equation [53] encountered in 
investigating scaling properties and general morphology of nonequilibrium systems. One possibility 
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of dealing with such a case is of course to compute directly the cumulants and construct the KM 
equation as described in section 2. A procedure similar to the one described in this section is 
technically more demanding and we leave it for future work. 



3.1 A useful check: the hnear case 

Before concluding this section, it might be interesting to assess the validity of the adiabatic procedure 
followed so far. To do so, we reduce ourselves to the linear case, i.e. when the potential V{V(^) is 
quadratic. Eq. ([37|) with the boundary condition ([39]) is exactly solvable 

PnniC, v^t) = i-^=L=exp (-l/2M^f - M^J^Cv^ - l^-l^v^) , (63) 

where 



M.,., = l^i^-e-n- (64) 

It is clear from this expression that the typical time scale by which the variable reaches the 
equilibrium value is the inverse of the Hubble time. At large times, t ^ H~^, the probability 
reduces to 

1 f27HA'^' 

PnniC,vc,t) = — ( ^ 1 [-m^/DctK^] ewm/D^t)Cv^]exp [-{6H/D^)vl] . (65) 

Integrating this expression in one finds that at large times 



and 



PUCt) = -^^_exp [-{9H^/2D^t)e] (66) 

{(') = I dCePun{C,t) = ^t. (67) 

Let us now check that the same result is obtained using the elimination of the fast variable V(^. In 

the linear case the equation for reduces to that of the harmonic oscillator and therefore A„ = 3Hn 
and 

^^'Hvd = 2^!^" (V^^/^C^c) exp {-3Hvl/D^) , 

^^^^Hvd = ym/D^v^ , (68) 
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where Hn{x) are the Hermite polynomials. We now proceed computing the operators Lo,0) -^o,™ 
and Lo,n- The operator Lo,o is clearly vanishing. To compute -Lo,n) we make use of the fact that 
y/3H/D^V(; = = 1/27^1 {y^3H/D^V(;) . Using the orthonormality condition (gS]), this gives 



This means that we need only to compute Li q, which can be done by exploiting the fact that 



il,o = -|;§. (70) 



QH/Dc_vl = 1/2'!'^°^^ + 1, giving 



Therefore the equation satisfied by the slow variable Q, from (j53p is, the KM equation 

^^'K.*) = I^^^K.') + -. (71) 

which gives (C^) = {D(^/9H^)t and coincides with (j67p as expected. Notice that the same result 
can be also obtained in an alternative way. From Eq. (j37p with quadratic potential V{v,^), we may 
derive the set of equations 



^^{vD = -m{vl) + D^. (72) 



At times larger than a few Hubble times, the correlators involving decay promptly to their 
equilibrium values 

resulting in 

-ie) = ^ (74) 



whose solution reproduces once more (j67p . 



4 Coupling to an en extra scalar field 

As further example, we consider the inflationary theory originally studied in Ref. ^ where an inflaton 
with a standard kinetic term is rolling down a flat potential and is interacting gravitationally with 
scalar fields cj„. We will consider from now on only one single scalar field as our conclusions do 
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not change if the number of fields is larger than one. We will also suppose that the scalar field has a 
potential V{a) and we will focus on its consequences on the probability for the comoving curvature 
perturbation as done in Ref. [58] and, partially, in Ref. |10]. The action we are concerned with to 
extract the IR behaviour therefore becomes 



S 



I 



1 + 



1 d 



(75) 



where the dots stand for further interaction between ^ and a that we disregard and the interaction 
between and the scalar potential V{a) acquires that particular form once we take into account 
that the linear dependence in the metric, ds'^ = (1 + C,/H)'^(lt^ — a?'{l + ZQdy? . Following again 
the rules dictated by [35] , the IR part of the comoving curvature perturbation and the field a satisfy 
the equations 



(76) 



a + 3H& + V'{a) = ry^ , 
where we have introduced the Gaussian and Markovian noise rja with correlator 

{r,^{^,t)vA^,t')) = ^6D{t-t') . (77) 

In Eq. (j76p we have neglected for simplicity the term with two derivatives with respect to time of 
the comoving curvature perturbation; furthermore, the (" dependence in the equation of a disappear 
at the order we are working upon the substitutions dictated by Ref. [35]. Introducing the variable 
Va = a, Eq. ([76|) can be rewritten as 



v^ + 3Hv^ + V'{a) 
The equation for the probability P{C,a,Va,t) then reads 



(78) 



d_ 
dt 



{3Hv^ + V'{a)) d_ 



da 



with initial condition 



8^2 dy^ 



(79) 



(80) 
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The operator 



can be written as 



L„. = -exp - 1.= !^) (3iJ„t6 + W + 6t6) e.p Q,^^ + , (82) 



where 



6t 



2V6TT dv„ 3i/2 



The normaUzed eigenfunctions ^n{va) of the operator 6^6, ie. 6^6^'n(f(T) = n^n{va) with > 0, 
are the eigenfunctions of the standard harmonic oscillator and can expressed in terms of the Hermite 
polynomials. They can be used to expand the probability as 

P(C, a, Va, t) = ^'oK)exp (-47rV(c7)/3iJ^) ^ P„(C, a, t)^n{va) , (84) 

n=0 

where 



^oK) = (-2^2^2/3^4^ 3i72/87r ^3^ . (85) 
The functions Pn{C,,(J,t) satisfy the equations 

- 6g2gM2 247r2 aC ( - l)^"-2 + + V(" + l)(" + 2)Pn+2 + (« + 1)P„) 

PI ^ 

^ 18i72 5C2 " 

- ^/^DPn-l - SHnPn - . (86) 
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They split into 



d _ _ V'{a) 3H^ d 3H d 



,2 



and 



V'(a) 3H^ d / ^ \ - 

where we have supposed that t S> so that the probabihties P„ with n > 1 have asymptotically 
reached stationary solutions. It is easy to convince oneself that the cr-dependence of Pq is of the 
stationary type 

Po(C, t) = Po(C, t)AAoexp (-47rV(a)/3i?4) , (89) 
where A/q is a normalization factor. The expression above satisfies DPq = 0. Due to the orthonormal- 
ity conditions for the ^{vcr), the probability P{Cit) coincides with Po(C)^) we get the following 
KM equation for P(C,i) 

8 1 /)2 1/93 

^PC.t) = 2°f8?^K,t) + 3fSfg^i'(C,t) + -, 

where one recognizes that the corrections to the tree-level diffusion coefficient originate from the 
one-loop contribution to the two-point correlator of Q obtained by exchanging the a through the 
potential V{a) (second term in the right hand side of of the middle line) and a = (last term in the 
right hand side of of the middle line) degrees of freedom. The term proportional to one derivative 
with respect to C, is absent once it is insured that the expectation value of C, vanishes at all orders. 
The fact that ((^^) ~ t (again neglecting the time change of the slow-roll variables) is due to the fast 
dynamics of the a field and which can be integrated out as they reach equilibrium values on a 
time scale of the Hubble time. 



5 The case of interactions dominated by slow-roll terms 

As last example, we will be concerned with the case in which inflation is driven by a single scalar 
field with canonical kinetic term. In such a case the self-interactions of the comoving curvature per- 
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turbation are dominated by terms which are suppressed by slow-roll parameters. To (considerably) 
simplify the matter we work in the limit in which the slow-roll parameter \rj\ = \e/eH\ S> |e|. This 
happens if the model of inflation is, for instance, an hybrid model with potential V{(j)) = Vb+m^</>^/2; 
in this case r] is constant and equal to {tvP/SHq) = (m^Mpj/Vo), while e = (j)fexp{—r]N) (where 
is the number of e-folds till the end of inflation and (pf is the value of the inflaton field at the wa- 
terfall transition ending infiation) is exponentially suppresed jl]. The curvature perturbation action 
simplifies to [57] 



S = M^i I (fx 



dt 



2dt 



91) 

Applying the rules of Ref. |35j and working to order 0{r}) we find the equation for the IR part of C 



C + 3HC-3Hr]CC = 3H7]^. 
With the usual position C = previous equation can be written as 



(92) 



V(; + 3HV(;{1 



C 



(93) 



At this stage we might proceed as in the previous section by integrating out the fast variable 
after having found the eigenfunctions and eigenvalues of the corresponding operator L^^. We take 
this time a short-cut and eliminate as follows. In Eq. (j93p we neglect the time derivative of 
and solve for the latter 



1-7?C 



which gives 



c 



(94) 



(95) 



l + r]C 

Instead of finding the corresponding KM equation, we follow this time an alternative road. We 
integrate directly Eq. (|95p after having changed the variable of integration from ttoN = f^ H{t')dt' 



C{N) = - 
V 



exp 7] 



N 



1 



(96) 



Suppose we are interested in the variance computed at some e-fold which corresponds to a scale 
inside our current comoving Hubble radius. Therefore we can write = (A^ — Nqq) + A'^go where now 
A'^go corresponds to the number of e-folds computed from the beginning of inflation till the moment 
when the wavelengh corresponding to our current Hubble radius exited the horizon during infiation. 
We can therefore write 
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C(iV) 



exp (rj [ dN'^+r] [ dN' ^ ) - 1 
\ Jo H J^^^ H J 

exp(,r"d'V'|)(l+.r/A"|)-l' 



(97) 



Due to the fact that the noise is white and therefore there is no cross-correlation between noises 
evaluated at different times, the variance becomes 



(C^, . (e.p (, f " diV' I) exp (, f " d^V" I) ) ( £ diV' | £ diV" |) , (98) 

2\ 



As the noise r/^ is Gaussian and for any Gaussian quantity O the property (expO) = exp 2(0 
holds, we obtain 



(C^(iV)) ^ exp I i!^l_il!(c2(iVgo)) ) (c^(iv - ATgo)) , (99) 



where we have traded the slow-roll parameter 77 with the spectral index of the comoving curvature 
perturbation [n^ — 1) = 2%]. This result coincides with the semi-classical result obtained in Ref. 

once one neglects there the contribution from the running of the spectral index which in our 
example is ^ . Indeed, in Ref. pT] the power spectrum for modes inside our current Hubble 
radius H^^ receives IR contributions once it is averaged by changing the location of the sphere with 
radius ~ Hq^ inside the superlarge box (corresponding to the total number of e-folds N). This is 
not surprising since the power spectrum computed in our Hubble volume, but averaged inside the 
superlarge box, is nothing else that the power spectrum inside the superlarge box. The latter is 
affected by large IR divergences, but is not directly observable. 

Notice that the result can be obtained by writing the corresponding approximated equation 
governing the probability P{C,t) |56] . The adiabatic elimination of the fast variable adopted in 
the previous sections indeed gives 



d „ Dt; d 



(1 + 7/0^(1 +r/C)P(C,iV) 



where dX = dC/(l + ?/C)- Recalling that P{X, N)dX = P(C, N)dC and to the order in the slow-roll 
parameters we are working at, we obtain that 

{C\N)) j dxl [exp {i^X) - if P{X, N) . (101) 

As P{X, N) is approximately the Gaussian function (after a change variable for N to incorporate 
the time dependence of the diffusion coefficient), one reproduces the result ([99]) . 
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6 Conclusions 



The study of interacting fields in cosmological backgrounds is a field which has attracted recently a 
lot of attention as future experiments will measure correlation functions of cosmological fluctuation 
variables with large accuracy. Cosmological correlators will be affected by ultraviolet divergences, 
as in quantum field theory, and by infrared divergences from modes which left the horizon at some 
time earlier than the relevent observable mode. The stochastic approach offers a technique to study 
such infrared effects. The infrared part of the curvature perturbation performs a random walk under 
the action of the stochastic noise generated by those modes which continuously leave the Hubble 
radius. The equation governing the probability distribution of the comoving curvature perturbation 
to acquire a given value at a given spatial point and time is the Kramers-Moyal equation, which 
generalizes the more standard (at least in cosmology) diffusive Fokker-Planck equation. In section 2 
and 3 we have for simplicity ignored slow-roll suppressed interactions, noting that the interactions 
leading to large IR effects are slow-roll suppressed and hence ignored in this limit. In section 4, 
we have shown that in the specific example, considering the interactions with a test scalar field, 
the late time behaviour of the effects on the two-point correlator of the curvature perturbation ^ is 
consistent with the fact that C, is conserved on superhorizon scales. The infrared effects considered 
in the example do not alter this property, as expected. In section 5, we considered a specific 
example of slow-roll inflation, and were able to reproduce the results of the semi-classical relations 
in [11]. Our work is far from being completed. First of all, we have derived in its full generality 
the Kramers-Moyal equation from a path-integral approach and shown how to alternatively derive 
such an equation through the adiabatic elimination of fast variables in a selected set of examples. 
One would certainly like to generalize and systematize this procedure. In particular, the role of self- 
interaction terms where spatial gradients of the curvature perturbation are present have not been 
addressed in this paper. This is because stochastic approach turns out to be surprisingly technically 
cumbersome when including derivative interactions. But although a simpler alternative semi-classical 
method has been developed in p4] (see also [271 [29] for further discussions) for computing IR 
corrections to correlation functions of cosmological perturbations, the stochastic approach may still 
be useful for understanding some alternative or complementary questions. Furthermore, we have 
taken for granted that the stochastic approach leads to a resummation of the leading infrared effects. 
While results in the literature based on the study of self-interacting scalar fields in a de Sitter stage 
support this expectation, a more formal derivation for the case of the curvature perturbation and 
its interactions is desirable. 
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